Force and flow transition in plowed granular media 
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We use plate drag to study the response of granular media to localized forcing as a function of 
volume fraction, <f>. A bifurcation in the force and flow occurs at the onset of dilatancy, <f> c . Below 
<j> c rapid fluctuations in the drag force, Fd, are observed. Above <f> c fluctuations in Fd are periodic 
and increase with <f>. Velocity field measurements indicate that the bifurcation in Fo results from 
the formation of stable shear bands above <j> c which are created and destroyed periodically during 
drag. A friction-based model captures the dynamics for cj> > cj> c . 

PACS numbers: 47.57.Gc, 45.70.Cc, 47.20.Ft, 83.80.Fg 



Granular materials fascinate because they can act like 
both fluids and solids Recent work has focused on 
the static problem of mechanical rigidity (jamming) in 
which the packing density <f> (the ratio of solid to occupied 
volume (2|) is increased until grains crowd sufficiently to 
develop a finite yield stress |3|]. Less work has explored 
the related process of "unjamming" Q where initially 
jammed granular ensembles flow in response to forcing 
and where the initial packing density plays an important 
role: varying <f> changes the local packing structure of 
grains which in turn affects the flow and force dynamics 
of the material response. 

We are interested in granular media subject to local- 
ized forcing (for instance from limbs @). In general, 
granular systems sheared at the boundaries evolve to a 
steady-state <f> c Q: depending on initial <f>, the medium 
compacts (4> < </> c ) or dilates (<j> > </> c ) as cj) — > </> c . In 
contrast, localized forcing, realized by an intruder trans- 
lating through an initially homogeneous medium 0, H| 
and viewed in the reference frame of the intruder, drives 
material only in the vicinity of the intruder toward <p c 
while simultaneously advecting undisturbed media into 
the flowing region. The result is the repeated unjamming 
of "fresh" material and the possibility of complex spatio- 
temporal dynamics in and around the zone of disturbed 
material surrounding the intruder. 

Here we drag a flat plate through initially homogeneous 
granular media at prepared (j) to continuously drive the 
system away from (fi c and measure the resulting force and 
flow fields. As <j> increases, the onset of shear dilatancy at 
4> c drives a bifurcation in force dynamics and media flow. 
This bifurcation is governed by a complex spatiotemporal 
flow response associated with the stability (instability) of 
shear bands nucleated by the plow above (below) <j) c and 
their evolution in response to plow motion. A model of 
shear band dynamics captures the oscillatory behavior 
above 4> c and suggests complex dynamics below <p c . 

Experiment- Drag measurements, see Fig. l(a,b), 
were performed in a 27x86x13cm 3 bed of polydisperse 



256±44/mi glass beads (Potters Industries; density p — 
2.51 gcm~ 3 ). Similar effects to those described here were 



observed in other granular materials (see EPAPS [12j), 
including heterogeneous beach sand and poppy seeds. 
Air flow through the porous floor initially fluidizcd the 
medium and then a combination of air flow (below flu- 
idization) and mechanical vibration generated the desired 
initial volume fraction (0.579< <j> <0.619). Air flow was 
turned off prior to testing and volume fraction deter- 
mined from bed height images as <j) — M/pAh, where M, 
A and h are the bed mass, area and height respectively. 
A stepper motor and linear translation stage displaced 
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FIG. 1. Experimental overview: (a) drag force, surface de- 
formation and (b) velocity fields are measured as a function 
of prepared volume fraction, <j>. (c) Mean drag force increases 
linearly with <j> while (d) temporal fluctuations in Fd become 
periodic as <f> is increased (arbitrary vertical shift). 
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a 3.9 cm wide by 0.3 cm thick steel plate submerged to 
a variable fixed depth 7.5 < d < 9.5 cm over a distance 
of 50 cm at a constant speed of v = 4cms~ 1 . An opti- 
cal encoder recorded the position and a 3-axis ATI load 
cell mounted between plate and translation stage mea- 
sured the drag forces (sampled at 200 Hz). Using laser 
line profilometry, we recorded the resultant surface pro- 
file and used it to quantify the change in vertical cross 
sectional area, AA, normalizing by the submerged area 
of the plow, Ap — width x depth. Profiles were mea- 
sured 35 cm from the start of drag where the profile was 
in a steady state. In separate drag experiments a flat 
plate at depth 5.5 < d < 6.5 cm was positioned against 
a transparent wall and displaced at v = 2 cms -1 par- 
allel to the wall to image grain flow. High speed video 
(250 fps) of the flow was recorded and analyzed in Mat- 
lab using image registration with a correlation time step 
of 0.02 s 0. The near- wall setup exhibited similar force 
fluctuations as the bulk but with approximately half the 
mean force. We remove two systematic variations in the 
force — a slow increase in the baseline force during drag 
(« 5% of the mean) and the decrease in depth of the 
constant height intruder with increasing <f> — by defin- 
ing the drag force, Fd, as the raw drag with slow drift 
removed (3rd order polynomial fit subtracted while pre- 
serving the mean) multiplied by a depth correction factor 
(^d^) 2 normalized to the loose pack depth dp p. Separate 
measurements at controlled intruder depth support this 
normalization technique. 

Force bifurcation- The mean drag force increases ap- 
proximately linearly with <fi as expected due to in- 
creased bed density and average coordination number, 
see Fig. 1(c). A bifurcation, however, occurs in the force 
fluctuations: Fd at lower <f> exhibits small amplitude, 
fast, and <f> independent fluctuations, while at higher </>, 
slower, periodic, and larger amplitude oscillations in Fp> 
occur which grow in duration and magnitude as <fi is in- 
creased, [Fig. 1(d)]. 

We characterize the bifurcation by measuring the stan- 
dard deviation of force, o~f d , and the average plow dis- 
placement between positive slope zero crossings of the 
force, AL, as a function of initial <f>. As <f> increases we 
observe a sharp bifurcation in aF D and AL [Fig. 2(a)] at 
the compaction/dilation transition (i.e. AA/Ap = at 
(j> c = 0.603 ± 0.003), see Fig. 2(b). As the bifurcation 
is approached from below, ap D is small and constant; 
above </> c and with increasing cj>, ap grows linearly while 
AL increases nearly linearly For (f> > <f> c , AL is indepen- 
dent of plow speed (2 < v < 8 cm/s), see Fig. 2(a) inset, 
revealing it as a characteristic spatial scale and implying 
that the temporal dynamics of Fd result from granular 
flow mechanisms that change with <f> across the dilation 
transition. 

Flow bifurcation- Direct evidence of a bifurcation in 
the flow is revealed by observations of the displacement 
of material at the surface near the plow [Fig. 3(a)]. Below 
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FIG. 2. Bifurcation in force occurs at the onset of dilation 
4>c = 0.603 ± 0.003. Small and nearly constant for low <f), (a) 
standard deviation of Fd grows linearly above 4> c while (inset) 
force oscillation length scale increases nearly linearly (2, 4, 6, 
8 cms in red, blue, green, and black), (b) Compacting and 
dilating response of plowed granular packings vs. initial <j>. 
(Inset) Trough cross sections at high (black) and low (red) <f>. 



4> c surface deformation is smooth, and the boundary be- 
tween upwardly moving grains and the undisturbed sur- 
face advances uniformly at a fixed distance ahead of the 
plow. At and above <p c , however, the surface flow takes 
the form of periodic radial upwellings of grains which 
give the region around the plow a stepped appearance. 
The generation of surface ripples is correlated with fluc- 
tuations in Fp> and indicates that the bifurcation in ap D 
and AL results from a change in grain flow dynamics. 

To gain insight into the nature of the bifurcation we 
compare the grain velocity field at low and high <p in a 
vertical plane perpendicular to the displacement of the 
plow, see Fig. 3(b) and [l0| . In both cases, flow is largely 
confined to a wedge-shaped region with angle 9 in which 
grains move forward and upward in advance of the plow. 
At both (f>'s a shear band separates the flowing wedge 
of grains and the effectively solidified grains outside the 
wedge; the lower end of the shear band terminates at 
the bottom of the plate while the upper end reaches the 
surface. 

Temporal flow dynamics- Despite these common fea- 
tures, differences in spatiotcmporal evolution exist be- 
tween flows at low and high (f>. At low flow is inter- 
mittent [Fig. 4(a)]: the horizontal extent of the flowing 



FIG. 3. Surface flow and velocity fields differ below <f> c (left) 
and above <f) c (right), (a) Surface deformation from a sub- 
merged plate shows step like features above (j> c . (b) Snapshots 
showing differences in velocity field structure from transpar- 
ent wall experiments [Fig. 1(b)] 



region, 77, (measured at approximately half the intruder 
depth, 3.2 cm) rapidly advances and retreats as the plow 
moves into new material. At high 0, the spatiotemporal 
evolution of the flowing region is periodic with 7/ remain- 
ing fixed for long periods of time until rapid, repeated 
fluctuations precede a jump forward to a new location. 
We quantify the pinning of this flow boundary as the 
fraction of time 77 remains stationary (-5 = 0) during 
the entire drag. As <f> increases rj becomes increasingly 
stationary indicating that the spatial shear regions are 
more stable. Above (j) ks <f> c , the shear band remains sta- 
tionary for more than 50% of the drag duration. 

Shear zone localization can be understood by the pro- 
cess of granular shear weakening/strengthening. Below 
</> c material compacts and strengthens under shear, con- 
tinually frustrating the development of fixed shear zones 
and forcing shear to occur along a constantly changing 
set of failure surfaces. Above 4> c however, material di- 
lates and weakens under shear, causing flow to localize 
along a fixed plane. In a plowed system, shear localiza- 
tion for initial 4> > 4>c holds only shortly after shear band 
formation: as the plow advances into undisturbed ma- 
terial the localized band is forced to adjust. Using the 
average angle of the flow boundary, 0, to characterize the 
orientation of the shear zone [see Fig. 3(b)], we find that 
below 4> c force and flow are largely uncorrelated while 
above tfi c , 9 and Fd are strongly correlated with Fd and 
8 increasing in concert, see examples in Fig. 5(a). 

Wedge flow model- We develop a simple model based 
on the correlation of Fd and 9 to gain insight into the 
dependence of Fd on the flow state and <j>. As the bot- 
tom of the shear zone is pinned to the lowest point of 
the plow, the flowing region is taken to be a triangu- 
lar wedge of mass m — 2tane sliding up a plane in- 
clined at angle 9. Kinetic friction between the wedge 
and plane, n((j>), is assumed to increase with <j>. Bal- 
ancing the forces from the plow (assumed horizontal), 
gravity, and sliding friction, the model predicts the plow- 
ing force F(0,/x) = , where W is the plow 
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FIG. 4. Temporal evolution of flow differs below and above <j> c . 
(left) Space-time plots of flow speed along a constant depth 
horizontal strip show that the leading edge of the flow, 77, (a) 
fluctuates in loose pack and (b) is stationary between shear 
band formation events (four shown) in close pack, (right) 
Time series of |2 highlight chan ges in stability, (c) Percentage 
of time during 5 s drag that the flowing region is stationary 

(|f = 0) increases with (j> and is fit to [l + e" (< ^'M . 



width [Fig. 5(c)]. The plowing force diverges at 9 = 
(infinite block) and 9 = tan -1 and is minimum 

at 9 m i n = tan -1 /i 2 + 1) — fij . 9 m i n decreases and 
F(9 m i n ) increases with increasing fi((p). 

The model predicts that an initially jammed and 
homogeneous material with volume fraction O shears 
at an angle 9 m i n when the force reaches F min ([i), see 
Fig. 5(b,c). For O > tf> c , shearing along the slip plane 
dilates and weakens the material locally, causing fi(4>o) to 
decrease to ^(<p c ) which reduces Fd and slightly increases 
9, see path A in Fig. 5(c). As the plow advances the an- 
gle of the weakened shear zone gradually increases along 
F[9, fi(<p c )], indicated by path B. The shear band remains 
fixed at the surface (causing 9 to increase) instead of ad- 
vancing with a constant angle because, evidently, less 
high <j) material needs to be broken in the former case. 
With increasing 9, the force to push the wedge 9 max even- 
tually equals the force required to break the stronger ma- 
terial at <j) > <j> c , i.e. F[6 max , /i(0 c )] = F[6 min , fi((p )]. At 
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this point a new shear zone forms in front of the plow, 
decreasing 9 to m in(4>o) and the cycle repeats (path C). 
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in the surrounding bulk. With the continuous advection 
of undisturbed weaker material and the heterogeneity in 
4> in the previously sheared region, predicting the shear 
process below <p c is more challenging. 

In summary we observe a flow and force bifurcation 
in localized granular drag which occurs at the onset of 
dilatancy. Localized forcing in granular media induces 
a heterogeneous <f> distribution in the region of flow and 
represents a new way to test theories of granular flow 
[la |. Previous studies of drag in granular media Q were 
performed at slow speed to observe stick-slip fluctuations 
over distances on the order of a grain diameter; the ef- 
fects of varying initial <fi were not investigated. Our ob- 
servations of a bifurcation in (f> support the findings of 
Schroter et al. fl3| who reported a transition in dynamics 
at <f> = 0.598, and our measured <j) c is similar to that re- 
ported (4> c = 0.591) in recent impact experiments It 
is known that 4> c is a function of the applied stresses and 
cohesive / frictional interactions between the grains 0, E3] 
so we do not place importance on the specific value here. 

We thank Harry Swinney and Mathias Schroter for 
stimulating discussion, and Andrei Savu for apparatus 
fabrication. This work was supported by the Burroughs 
Wellcome Fund and the ARL MAST CTA under coop- 
erative agreement number W911NF-08-2-0004. 
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FIG. 5. Force dynamics and wedge model of granular flow, 
(a) Wedge flow angle 6 and Fd are correlated above <j> c . (b) 
Sketch of wedge flow model dynamics above 4> c with arrows in- 
dicating motion of shear band for the three phases A,B,C also 
shown in (a,c,e). (c) Dynamics predicted by model (arrows) 
showing family of (6,Fd) curves for varying n(cj>) and tran- 
sitions between fJ,(cj>) curves corresponding to 4>o (red curve) 
and <j)c (blue curve), see text. Comparison of model (9,Fd) 
curves and experiment: (d) for tf> < <f> c experimental data 
are bounded by model curves of arbitrary /j, < fJ,(4> c ); (e) for 
4> > 4> c model shows good agreement with experiment. 

The dynamics of the experimental data for (j> > <j) c 
are captured by the model as seen in (8,F(9)) space 
[Fig. 5(e)] and Here the system evolves between 

two curves [fi = 0.65 (blue) and [i — 0.78 (red)] and 
the transitions between these states qualitatively match 
the predictions of our model. The model predicts the in- 
crease in force fluctuations with increasing <fi above <p c . 
However for (f> < c/) c the evolution of Fd with 9 is more 
complicated because of the lack of stability of the sheared 
regions through shear strengthening. As we assume a 
feedback process exists between the flow dynamics and 
Fd, strengthening of the material during shear results 
in the relocation of the shear band to a weaker region 
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